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Abstract. The processes of the single Vy-production in hadron-hadron collisions are 
suggested for investigation of the nucleon spin. An approach is proposed for the de- 
termination of quark spin densities at low x. The lowest order electroweak radiative 
corrections to the observable quantities are calculated. The numerical calculations of 
the cross sections and the single spin asymmetries taking into consideration the elec- 
troweak corrections at RHIC energies have been made. 
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O . 1 Introduction 

^^ , Despite remarkable success in solving the proton-spin problem, which was risen 

f| ' by spin crisis [1], in experiments with fixed targets at CERN [2], SLAC [3] and 

Qh, HERMES [4] a conclusive solution to the problem of determination of u— , d— , s— 

Q.' quark and gluon polarizations has yet to be found. Purely inclusive measure- 

D , ments determining the spin structure functions for nucleons and deuteron are 

i-»5 ' unfortunately restricted to probe only certain combinations of the polarized par- 

ton densities. A full analysis would require additional inputs from other measure- 
ments to separate different components. For example, semi-inclusive measure- 
^ , ments with SMC [5] and HERMES [6] allow to access to a variety of polarized 

j^ ' parton distributions. However, direct and separate measurements of the polar- 

ized parton distributions will remain a limited significance for some time. 

Some new combinations of the polarized quark densities can be obtained from 
the experimental data on single spin asymmetries, when the polarized target is 
used in conjunction with an (un)polarized proton beam. This opportunity is 
offered by the dedicated RHIC spin program at BNL [7] which is supposed 
to start in a short period of time. Collider polarized experiments STAR and 
PHENIX at the RHIC energy {^/S « bOOGeV) wiU be the point where a high 
energy nucleon-nucleon spin physics will be studied by measuring a variety of 
spin asymmetries. Besides, there are some prospects on an acceleration of the 
polarized protons to ITeV at Fermilab Tevatron [8] and these will lead to similar 

physics, which will be accessible to RHIC and possibly to HERA-iV [9]. 

In this paper we are concerned with experiments which may provide direct 
measurements of new independent combinations of the quark densities in polar- 
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ized nuclcon. We will focus on inclusive single W^-boson production in hadron- 
hadron interactions with one longitudinally polarized beam 

A^i + N2-^ W^ + X -. l^ + X. (1) 

It's well known that the single spin asymmetry of the production of T4^-bosons 
via PP^> W^X (e.g. Al from Ref.[10]) is sensitive to the form of polarized quark 
distributions. In Sec. 2 of our paper we get formulas for single asymmetries 
of the process (1) for the case when in the final state only charged lepton is 
detected. Moreover, we will show the possibility of proposing reactions to study 
the polarized quark densities in the region of small x. 

Radiative events originating from the loop diagrams and (as we consider 
the inclusive process) the processes with the emission of real photons, cannot be 
removed by experimental methods and so they have to be calculated theoretically 
and subtracted from measured cross sections of proposed reactions. In this paper 
the total 0{a) electroweak radiative corrections (EWC) to the cross sections 
and single spin asymmetries have been calculated. In Sec. 3 and 4 we present 
the contributions of additional virtual and real emitted particles respectively. 
Numerical analysis and conclusions can be found in Sec. 5. Some technical details 
of calculations are presented in Appendices. 

2 Born cross section and single asymmetries 

The cross section for the inclusive hadronic reaction (see Fig. 1,1) is given in 
Quark Parton Model (QPM) by conventional formula 

'^%,^.^.±x=/^"i^"^ ^ /«''-^(xi,Q^)/f^^^(x2,02)rf4, (2) 

i,i';rl,r2 

where /» (a;, Q^) is the probability of finding constituent i with the fraction 
X of the hadron's momentum and helicity r (at squared momentum transfer in 
the partonic reaction Q^) in hadron a and da^w is the cross section for the 
elementary process leading to the final state (see Fig.l). The summation runs 
over all contributing parton configurations and over helicities of first and second 
partons (ri^2 = ±1)- The meaning of operator "hat" will be explained below. 

The process of the single VF-boson production in hadron-hadron collision can 
be described with a very good approximation by two pair of quark- antiquark 
subprocesses. So for the W^ production we have 



q2{pi) + qi'{p2,V2) - 


-^w-{q)- 


-.l-{h) + 9{k2), 


(3) 


qi{pi) + qt'{p2-,m) - 


^w-{q)- 


-.r(fci) + p(fc2). 


(4) 


and for VK+ one 








qt{pi) + qi'{p2,V2) - 


^ W+{q) - 


-.l+{k,)+l^{k2), 


(5) 


qi{pi) + qi'{p2,V2) - 


^ W+{q) - 


^l+{k^)+v{k2). 


(6) 
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Our notations arc the following (sec Fig. 1,2): pi is a 4-monrentunr of a first 
unpolarized (anti)quark with flavor i and mass mi; p2 is a 4-momentum of a 
second (anti)quark with flavor i' (i' denotes the weak isospin partner of the 
quark i), mass 7712 and polarization vector 772; fci is a 4- momentum of a final 
charged lepton l^ or /+ with mass mi; fe is a 4-momentum of (anti)neutrino, 
q = Pi + P2 is a 4-momentum of the H^-boson with mass mw . We use the 
standard set of Mandelstam invariants for the partonic elastic scattering 

S=(Pl+P2)^ t=(pi-fci)2, U^{ki-p2f. (7) 

Squaring the matrix elements of partonic subprocess and (we use the co- 
variant expression for the polarization vector of i' parton from [11]) we get the 
invariant parton-parton cross section in the Breit-Wigner form 



m,sis{{s-ml.)^ + ml.r^y^P'^P' "' ''^>ki,2k2, 



^4' - . .r „4 „//„ ^2\2 ,"2 r2 A p^ +P2-ki- fca) ,^ ^,^ , (8) 



^ii' — \ ,2/ 



where 

1/Nc ~ 1/3 is a color factor, s^ = -^/l — c^ is a sine of the weak mixing angle, 
Cw = m,w I'niz, Tnz is the Z-boson mass, Fw is the VF-boson width, Vn' is 
Cabibbo-Kobayashi-Maskawa mixing matrix, 

v?(\ — r2PN2) for (3) subprocess, 

t^(l — r2PN2) for (4) subprocess, 

t^(l + r2PN2) for (5) subprocess, 

v?{l + r2PN2) for (6) subprocess, 

and PN2 is the degree of longitudinal polarization of a second hadron such that 

PN2 = ±1- 

The integration w.r.t. 4-momenta of unobservable (anti) neutrino gives 

(fko 

— — S(pi + p2 - ki — k2) — 5(s + t + u - mf - ml — mf). (9) 

2K20 

Then in accordance with QPM we substitute Pi(2) ~^ a;i(2)-Pi(2): where Pi(2) 
is 4-momenta of initial nuclcons with masses tojv, a;i(2) is the fraction of the 
first (second) nuclcon momentum that is carried by the incoming quarks. We 
shall denote this procedure by operator " hat" . Then we multiply by the parton 
densities of first and second hadrons, sum over hclicities of quarks and integrate 
over xi, X2 according to formula (2). 

Let us introduce the Mandelstam variables for hadronic reaction 

S = 2P1P2, T = -2Pifci, U = -2P2fci, (10) 

(then s = X1X2S + mf + rrig, t = xiT + rairajq + mf , u ~ X2U + m2mN + mf) 
and integrate w.r.t. X2 with help of (5-function and noting that in QPM 

6{s + i+u-ml-ml-mf)^^S{x2 + '''^^^"'J'^^'^h , D = x,S + U + m%. 

We can see that in this case X2 = x^ = —{xi{T + mj^) + rrif)/D and this 
substitution corresponds to Born kinematics and we denote this by subscript (or 
superscript) "0". 
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Finally, let us consider the general form of the cross section of hadronic 
process (1). In the hadron-hadron collisions the center of parton-parton masses 
frame has an undetermined motion along the beam direction. Therefore we use 
the standard in this case variables: centre-of-mass energy (vo), component of 
the detected particle vector transverse to the beam direction (|fcij^| = fciT)j 
and pseudorapidity (77), which approximately equals rapidity y, since in our case 
Tni <^ kix] then we have for T and U 

T = -^/SklTe-\ U = -^/SklTe'' . (11) 

Integrating w.r.t. azimuth (p (it's possible as the first initial hadron is unpolarized 
and the second one is longitudinally polarized) we have phase space d^ki/kiQ => 
irdrjkirp, and hence the Born cross section has the form 

ao±=^ /"dxig,(xi,Q2)i:(,, (12) 

where a^ denotes the double differential cross section 

da -^ 
drjdkirp 

We use this abbreviation as well as the following 

a^ =a^ +PN2^cr^, (14) 

where a {Aa ) is the spin averaged (polarization) part of cross section in the 
whole text. 

The Born cross section is proportional to factor Sq = S{x2), 
where 

In the expressions (12), (15) the combination of the quark densities for a second 
nucleon has the form 

fP{x2,Q^) = qHx2,Q^) ~ c,,pN,AqAx2,Q^), (16) 

where {A)qi{x,Q'^) = f^{x,Q^){—) + J^{x,Q'^) arc (longitudinally polarized) 
spin averaged quark densities, Ci' = —1(4-1) for quark(antiquark). 

The form of the Born cross section (12) is convenient as the factorized part in 
the so-called 'soft' 0{a^) cross section (see Sect. 3, 4) and here let us present more 
simple expression for the Born cross section of the processes (1) as a convolution 
of partonic cross section with two parton distribution functions which we will 
use in the rest of this section 

(17) 
http : //link . springer . de/link/service/ j ournals/10105/index . html 



Electroweak corrections to the observables of H^-boson production at RHIC 



Let us define the polarization single asymmetries as normalized difference of 
the cross sections with specific spin configurations 

A (V, klT) = ^ TT- — +7 7T = -r^- (18) 

<^^{PN2 = 1) + cr±(pw2 = -1) cr± 

Supposing that CKM matrix has diagonal form Va' = Sw we can see that flavors 
of initial quarks and antiquarks are i = d,s,b; i' ~ u,c,t for (3) subprocess, 
i = u,c, t; i' — d,s,b for (4) subprocess, i = u,c, t; i' — d,s,b for (5) subprocess, 
i = d,s,b; i' = u,c,t for (6) subprocess. Neglecting the contributions of the 
heavy quarks (c, 6, t) (the strange quark contribution su —^ W~ is neglected too 
as we suppose Vsu equals zero) we can write the Born single asymmetries as 

;+ ^ Jdxiiu'ixi)Ad{x<i) ~ d'{x^)Au{x^,)) 

" ^''' '^^ /(ixi(u'(xi)d>0) + rf'(xi)u(xO)) ' ^ ' 

A{'n,kiT)^A'o{'n,kiT)iu^d}. (20) 

Here 

u'{xi) = Ktu{xi), u'{xi) ^ Ktu{xi), d'{xi) ^ Kud{xi), d'{xi) = Kud{xi), 

(we suppress the argument Q^ in the quark distributions of formulas (19,21)). 
The physically allowed region of xi and X2 (see Fig 2.) is given by 

-^^y^<.t<l, .^<.2<1. (22) 

Let us remark that in the region of large xi and small kirp/yS the expression Xj 
does almoust not depend on xi {x2 ~ —T/S). Dividing the region of integration 
in (19,20) by parameter x\ (we can choose such a value of x\ in order to have well 
defined polarized quark densities in the region xi < x\) we obtain the expression 

Au{-T/S) I dxid'{xi) - Ad{-T/S) f dxiu{xi) = (23) 



-4o"^('7,^it) / dxi{u'{xi)d{x'^) + d'{xi)u{x'^)) 



dxi{Au{x9,)d'{xi) - Z\d(a;^)u'(a;i)) 

and analogously for A'' by replacing u ^^ d. 

Off-shell approach, which we used gives the possibility to investigate the 
polarized quark densities in the region of very small x. Really, the expression 
—T/S under conditions of collider experiment can reach very small values, so 
using, for example, RHIC kinematics point \fS = 500GeV, kix ~ lOGel^, rj = 2 
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we can see that —T/S does not exceed 0.0027. That minimum value of x is 
about one order of magnitude below the values quoted in studies where on-shell 
approximation for the VK-boson was used. 

So, equations (23) connect the polarized quark densities in the region of small 
X with the observable single asymmetries, combination of unpolarized quark den- 
sities and polarized quark densities in the region where they are well defined. If 
the three of the supplementary measurable quantities (for example double asym- 
metries A±{x,y) from Ref.[12] and QPM-expression for .91(2;)) are used, equa- 
tions (23) allow to determine the low- a; behavior of polarized u— and d— quark 
and antiquark distributions in nucleon. 

3 Contribution of additional virtual particles to the cross 
section 

To extract the reliable data on single asymmetries with high precision from 
hadron collider experiment, it is necessary to consider higher order electroweak 
radiative corrections. The final state photonic corrections to processes of W- 
production in unpolarized pp-coUisions were calculated in Ref. [13]. More accu- 
rate calculation of these electroweak corrections have been suggested in Ref. [14] , 
where both initial and final state radiation have been included, and the coUinear 
singularity associated with final state photon radiation is regularized by the lep- 
ton mass. On the contrary, the coUinear singularities associated with initial state 
radiation are subtracted and adsorbed into the parton distribution functions. 
Due to these results do not contain the quark masses, but depend analytically 
on other unphysical parameter dg, which determines a coUinear region. 

In this paper we present the new explicit formulae for EWC to inclusive 
single M^-production in polarized hadron-hadron collisions, where by analogy 
to the calculations of the radiative corrections to the hadron current in the 
deep inelastic scattering of lepton by nucleons [15, 16, 17, 18] we used a finite 
quark masses to regularise the coUinear singularities. In Sect. 5 we consider the 
dependence of the results on various quark masses and the comparison of them 
with the results for EWC, when the strategy of coUinear singularity extraction 
has been employed. 

As in final state only charged lepton is detected we use for calculation the co- 
variant method [19], which has conclusive advantage: the results are independent 
of unphysical parameter - maximum soft photon energy Ecut (see for example 
Ref. [14]). 

The one-loop contribution of additional virtual particles (V-contribution) 
has been calculated in t'Hooft-Feynman gauge and in on-mass renormalization 
scheme which uses a, tovkj^t^z, Higgs boson mass mH and the fermion masses 
as independent parameters. The virtual one-loop diagrams are shown in Fig. 3. 

The cross section of V-contribution is proportional to the Born cross section 
and could be written as 

^v =E/'^^i'?«(^1'Q')^v U.=.o^o. (24) 
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Here factor Sy is polarization independent and consists of seven terms 

S't}' =Sw + Svi + 5vq + 5 SI + 6sq + 5^w + ^^w ■ (25) 

We do not repeat here the explicit expressions for all these contributions but 
instead refer the readers to Ref.[17],[20]. Let us consider the meaning of vari- 
ous terms in (25). The correction Sw is the VF-boson self-energy contribution 
(Fig. 3,1), which has a slightly modified form from what is presented in Ref.[17]. 
This is due to resonant W boson production 

c OD s - m'^ - imw Tw avk/ \ /oc\ 

^--2^^(7^<o^t^<;t^^-^^)' ^''^ 

where IJ)p (s) is the renormalized transverse part of the VK-boson self energy 
(it can be found from formula (6.1) of Ref. [20]). The term Svi is the leptonic 
vertex correction (Fig. 3, 2) (see formula (2.10) of Rcf.[17]), Sy is the quark 
vertex correction (Fig. 3, 3) (see formula (2.11) of Ref.[17], we left the flavour 
indices at this correction and boxes contributions) , Ssi is the neutrino self energy 
contribution (Fig. 3, 4) (see formula (2.12) of Rcf.[17]). Up-quarks self energy 
contribution (Fig. 3, 5) (down-quarks self energy contribution equals zero) is given 

by 






«Uln^--n^)-g^(ln^-21n^)4 



(27) 



This expression is obtained from formula (5.46) of Ref. [20] (see also the remark 
below it). 

We recalculated the box contributions through the quantities I12 and /f^ 
of Ref. [17] and present them here to stress their dependence on flavour of quarks 
and the channel of reaction: jW box contribution (Fig. 3, 6 and Fig. 3, 7) 

S:ii = S^tv^S-,w, S%^S%^S^w{t^u,t^i'), (28) 

where 

S^^ = —ReQi{QJ^'^'\s,t)+Q,,I^'^^''{s,u)); (29) 

TT 

ZW box contribution (Fig. 3, 6 and Fig. 3, 7) 

Sziv = ^zw = ^zw , <5|V = 4w ^Szw{t^u,t^t'), (30) 

where 

Szw = — i?e[((i;; +a/)(i;j + flj) + (wy +ay)(t;j/ + aj/))/f^(s,i) 

TT 

+ {{vi + ai){v,, + a,,) + {v^ + a^){v, + a,))/f '^(s, u)]. (31) 

The vector and axial vector couplings constants Vi, ai of the i-fermion~.Z 
vertex are determined from the electric charge (Qi) of fermion expressed in 
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tinits of proton's charge (e.g. Qu — Qu ~ +2/3, Q; = — 1) and the 3-coniponcnt 
of the weak fermion isospin [If) 



^i ^'^wWi ^i 



Let us present the cross section of V-contribution as the sum of infrared 
divergent (IR) and IR-finite parts 

4 = <^v'" + 4'^ = E /^^i*(^i' Q')^o('^(/' + ^f )U,=.o. (32) 

i.i' 

Infrared divergent part of V-contribution is regularized with the help of photon 
mass A and for the correction 6\^ we find 

'5y^-£log^J(0), (33) 

where the expression for J(0) will be considered in the next section. The finite 
part of V-contribution cross section contains correction 

5^ = 5X^ -6{,^ = 5^{\^^s). (34) 



4 The photon bremsstrahlung N-^ N^^^ l^^X 

In order to get infrared finite results for the hadronic process cross section we 
have to include the real bremsstrahlung correction (Fig. 4). 

The differential cross section for the partonic process with emission of one 
real photon reads 

pot 

where squared matrix element^ is the sum of initial state (index «), final state 
(index /) and interference terms respectively 

Y, \R? = (ARtT"' + (-Rz-R/ )""' + (AR^ + RiKT"', (36) 

pol 






^ For the calculation of matrix clement we used the standard Feynman rules and the proce- 
dure of separation of diagram Fig. 4, 4 into the initial and final state radiation parts according 
to Ref. [13]. As a result, 'initial' part of our matrix element completely coincides with formula 
(2.14) of Ref. [13]. 
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Indices n, n' denote the type of the processes: n — —(+) for qq{qq) initial state, 
and n' = —{+) for Z~j>(Z+j/) final state, and IIi{nq) is a T4^-boson propagator in 
the case of the lepton(quark) bremsstrahlung 

III = l/is-m'^ + imwrw), Hq = -l/{s-2kq-m^ + imwrw) = 1/iz + twr), 

(37) 
twr ^t^-imwRw, t^, ^ v - a + m^,q ^ Pi+ P2- 

As the kinematical variables of the radiated process we use in this case 

z — 2kki, zi — 2kpi, ti — {p2 — k2) , ui ~ 2kp2 ^ v + z — zi, 

V — 2kk2 — s + u + t — nil ^ ^^2 ^ '^; i (38) 

where fc is a 4-nionientuni of the radiated photon. 

The matrices U originate from the products of bispinor amplitudes and the 
expressions (1 ± 75): 

U^^p = (l-75)pi, Ur^ = {l+7r,m){p2-m2), !/+, = (1-75)P2, U^^ =pi~mi, 

(39) 
Ui,a = (1 - 75)^2, Ul^ ^ki+m, C/+„ = (1 - 75)fci, C/+P = fc2, (P - 7>m) 
and the matrices G originate from the fermion propagator and iyVF7-vertex: 

iG-r - Q.r'-^^^^^ - o.M^^,. + (crr: (40) 

(Gj ) ^Qi 7 -W7 1- (Gj ) , 

Zl Ul 

{Gjr = g.^^^7- + {Gf-r, (G+)- = -Qa-^^^ + {Grr, 

(GrU-^'''iO(-fc'?'fc-'?)' iGY')^^^l'''^C',.,,{-k,k-q,q), 



'MP 



2kq P'^P ' '"' -'' ^ / Vp ' 2fcg 



where the matrix C'^ corresponds to VFVF7- vertex and for the minimal standard 
model reads 

C^^p(p°,P+,P") ^ 9pv{p^ -P'^)tj,+9tj,v{p'^ -p")p+giip{p" -p^)v 

After transition to hadron-hadron cross section according to prescription of 
QPM, we can present the cross section of bremsstrahlung (R-contribution) by 
splitting it into a soft IR-part and a hard contribution [19] (we again use the 
abbreviation (13) for the double cross section a^) 



^R ^^R +^R ■ 



(41) 
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The infrared divergent part (the first term) of expression (41) reads (wc made 
the substitution dx2 = dv/D) 

'^fl''' = E/^^i*(2;i,Q2)^^fi, (42) 

where 

r^« = _- /"""^ dvS{x2)I[F''^]. (43) 

The procedure of integration over the photon phase space defined as I[A\ is 
described by (1) in Appendix A. And for F^^, using the method of separation 
of IR-terms [19] we find 

F^^^Q^^'4+ciQiQ~-ciQiQ,— + Q^^'4-Q,Q,^ + Ql^, (44) 

Z^ ZZl ZUl Zf ZlUl uf 



where 



_ , f 1 for (3) and (6) subprocesses, , , 

' —1 for (4) and (5) subprocesses. ^ ' 



Introducing J{v) = v hm /[F^^] we obtain two parts of the soft cross section: 



A^O 



Vmax Vmax 

r^« = --SoJiO) f — + - dv{SoJ(0)/v - ^(X2)/[F^^]). (46) 

n J V TT J 

The second term is infrared free and lower limit of integration for it equals zero, 
and the first one contains infrared divergence. In the center-of-mass-system of 
initial hadrons the limits of integration are 

^ D'^+TU 

Vmin = (fc= 0) = 2Afc2o = Au, V fH -=— , 

D\/ S 

Vmax = {X2max = 1) = D{1 - X^^) (47) 

and the infrared divergent part of bremsstrahlung cross section has the form 

-o/ttZ'o J(0) \0g{Vmax/'>^v). 

Summing up IR-parts of the V- and R- contributions (formulas (32) and 
(42)) we get 



cr^'^^ +(Ty'^^ "" E / '^^i'?«(^i' '9^)^^o^(0) log 



v'^so 



-j: fdxM-uQ')- 1 dv^^^^M^^Mm^ (48) 



b 



as a result infrared divergences completely cancel out. 
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During the calculation of I[F^^] the following expressions are used 

z^ mfv zj rriiV 

1 1 t^ 1 1 u2 

l[ — ] = — rlog^r^. I[ — ] = log^r^: 

zzi vt mfm^ zui vu mfm2 

I[^] = - log ^, I[^] = -^, (49) 

ziui vs m\m2 uj ttIjW 

and therefore J(w) has the form (we remind that v = D{x2 — x'2)) 

^2 ^ ^ , t 



J{v) = Qi - ciQiQilog^ — 2 +c;'9/Qj'log^ — 2 

mfm'i mfm2 

+ Ql - Q^Q^' log ^-^ + Ql . (50) 

?7lf?7l2 

After extraction of infrared singularity from partonic process cross section 
and integrating over whole phase space of real photon (see Appendix A) we can 

present the hard contribution cross section of inclusive process A'^i A^2^ / 7X 
in the form 

<yR^ = ^ f dXldx2q^ixl,Q^)F^^\x2,Q^)S^, (51) 



where 



a^ 



+Qinin+Re[Viu,] + Re[ni]Vg^ + nin+v^), (52) 



and 



77,77+ == l/((s - ml^y + ml^r^). 

The subscript in V means the origin of bremsstrahlung photon: l/q/w refers 
to the radiation from lepton (Fig. 4, 3) / quark (Fig. 4,1 and Fig. 4, 2) / VK-boson 
(Fig. 4, 4) legs respectively. Double subscript corresponds to the same interference 
term. The expressions for the quantities V can be found in Appendix B. 

So, we get three terms of 0{ct') cross section which include radiative correc- 
tions: infrared finite part of V(and R)-contribution cr^^o-, and the part [ay' + 
cr^' ), where infrared divergence cancelled out in the sum 

'^EWC="^V +'^R +\'^V +^R )■ (53) 

All of these terms depend on a second hadron polarization. This dependence is 
factorized in the expression 7^^, {x2, Q^), which is containing only in £'(2:2) and 
Eq. Such a factorization gives the possibility to obtain the cross section of a 
single ly-production in unpolarized nucleon-nucleon collisions by doing a simple 
replacment 

F^^P{x2,Q^)^qAx2,Q'). (54) 
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5 Discussion of numerical results and conclusions 

In the following the scale of radiative corrections and their effect on the observ- 
ables of the processes (1) will be discussed ^ . First of all we compare the numer- 
ical estimations of our formulas to the results of the calculations Ref. [14] for the 
unpolarized pp — > IX. We have calculated (Fig. 5) fi^ transverse momentum 
spectrum as a function of kij- for the process pp -^ ji^X at ^/S = l.STeV (Teva- 
tron kinematics). The integration over pseudorapidity (—1.2 < rj < 1.2) had been 
done and the MRS LO 98 set of unpolarized parton distribution functions [21] 
are used. It can be seen from comparison of Fig. 5 and Fig.7,b of Ref. [14] the 
total EWC have a similar behavior: 0{a^) cross section in the region 25GeV < 
kix < SAGeV is a little bit higher than Born cross section; the EWC contribu- 
tion to the cross section is insignificant at kix ^ 35(37)Gey for Ref. [14] (our) 
calculations; and in the region kirp > 39GeV the total EWC reduce the Born 
cross section in both approaches, so in the resonance region (fci^ ~ mvi//2) the 
difference between Born and 0{a^) cross sections is ^ 0. 004(0. 003)n6/G'eF for 
Ref. [14] (our) calculations. 

The main uncertainty for hadronic part of the EWC in our approach is the 
quark mass dependence. From Fig. 6 it can be seen that in the region kix < 
AOGeV the difference between corrections at various values of quark masses is 
rather essential: so increasing to„, md by a factor of 10 the correction decreases 
for ~ 0.05. In the vicinity of kij- ^ tovi//2 and for ki^ > m^/2 the discussed 
difference does not exceed 0.02. For the calculation of 0{a^) cross section of 
Fig. 5 and also for the rest analysis we used the values of the current quark 
masses of ?ti„ — 5MeV, md — 8MeV, considering the arguments concerning the 
choice of these not well-defined parameters, which are discussed in Ref. [16] (page 
1101, and references therein). However for quark masses of such an order a quite 
good correspondence with the results of calculations by the collincar singularity 
extraction method is reached (parameter 6g > 10~^). 

To estimate the scale of radiative effects and it's influence on observable 
quantities in the processes (1), the numerical calculations of the cross sections 
(Fig. 7 and Fig. 9), the total EWC S^^ to spin-averaged part of cross section 
(Fig.8) 

'^EWC = '^0* (1 + Sr± ) + PN2^4 (1 + S^± ) (55) 

and the single spin asymmetries A' taking into consideration EWC (Fig. 10) at 
typical values for future experiment STAR at the collider RHIC {^/S — 500 GeV, 
— 1 < rj < 2, A<P — 2tt) have been made. We used the following standard set 
of electroweak parameters: a = 1/137.036, mw = 80.43Gey, mz = 91.19Gey, 
Mh = 300.0Gey, the fermion masses: niu — 5MeV, nid = 8MeV, nis = 
IbOMeV, ruc ^ l.5GeV, wn, = 4.5Gey, m* = SO.OGeF; the GRV94 proton 
parametrization of parton distribution functions [22] for unpolarized quarks and 
the GRSV96 LO proton parametrization of parton spin densities [23] . We have 
selected for the Q^ in these Q^-dependent distributions as well as in Ref. [14] 

^A program of numerical calculations of the cross sections and single spin asymmetries of 
processes (1) (in FORTRAN) including 0(a) electroweak corrections is available by contacting 
to author via e-mail 
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Q = rn^, since trying to use for Q an explicit momentum transfer in the par- 
tonic reaction (e.g. Q^ = s for Born and final state radiation cross sections) 
we have not obtained some noticeable difference for the numerical estimation of 
the observable quantities. And at last, as the contributions to the observables of 
investigated processes which have different origin are not distinguishable exper- 
imentally, we did not conduct the comparative analysis for them. 

We present in Fig. 7 the numerical results for the spin-averaged part of the 
double differential cross section as a function of ki^p for the different values of 
pseudorapidity. As we can see, for e^ and e~ in the final state EWC are essential 
and increase the Born cross section in the region kij' < mw/'2 for r) = —1, 0; the 
EWC are not significant at small kirp for ?? > 1, and in the region kirp > m\Y /2 
for all values rj. For /i+ and /x~ in the final state EWC are most essential in the 
resonance region for rj ~ —1, 0, and at the small fciy in the case t] = 1.5, 2 
(EWC decrease the Born cross section). 

It can be seen from Fig. 8, where the corrections (5" to spin-averaged part 
of Born cross sections are represented, that muon corrections are lower than 
that of electron (difference is approximately 0.12). This fact corresponds to a 
common character of EWC dependence on masses which regularize the coUinear 
singularity (the mass rn^ is more than m^. by 2 x 10^ times) and which were 
discussed in the beginning of this section. General features of the corrections 
behavior for all cases (e , /i*) and for different t] are the increase in the region 
fci2- < m^^ 12, the sharp falling in the vicinity of resonance and smooth growth 
at k\j' > mw/2. The corrections to transverse mass distribution, which are 
represented on Fig. 9 of Ref.[14] have a similar behavior. 

The scale and the behavior of the radiative corrections to polarization part 
of the Born cross section are shown in Fig. 9, but as this part of the cross section 
is unobservable we shall not analyze it in details, and consider the influence the 
radiative effects to observable quantities - polarization asymmetries (19), (20). 
So, Fig. 10 shows Born asymmetries Aq together with the EWC corrected results. 
One can see that the corrections to asymmetries are significant and decrease the 
Born value practically in the whole investigated region except the region of small 
kix and large ry. Let us remark that the corrections are almost independent of 
particle mass in the final state {me or m^). In other words the ratios of the cross 
sections (i.e. asymmetries) are much less sensible to the values of masses than the 
cross sections themselves. For the same reason the corrections to asymmetries 
do not depend practically on the choice of quark masses too. So, if we change 
the masses of the light quarks from to„ = 5MeV, rud = SMeV to the extreme 
value rUu = md = O.SSGeF we get a shift of radiative corrected asymmetries by 
less than 0.01 in a whole investigated region, and in the region 77 = 2, kix = 
lOGeV, which is very important for the quark density analysis at low x (see 
Section 1) this shift equals 0.0032(0.0021) for ij,~{ij,'^) case. Thus we can suppose 
that our approach to the calculation of radiative corrections to the single spin 
asymmetries is practically free from the quark masses uncertainty. 

In conclusion we have presented the scheme to use the processes of the single 
W-boson production possible at hadron-hadron colliders for the investigation of 
the nucleon spin structure. An approach is proposed for the determination u— 
and d~ quark (and antiquark) spin densities at low x by a set of the observable 
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quantities. The electroweak radiative corrections to the differential cross sections 
and the single spin asymmetries of the processes (1) are calculated. The obtained 
formulas do not depend on any parameters of the infrared divergence extraction 
(e.g. maximum soft photon energy). We checked that EWC to asymmetries do 
practically not depend on the masses of the quarks which regulate in our ap- 
proach the coUinear singularity. Our results can be used for studies of nucleon 
spin at hadron-hadron colliders (RHIC and HERA) . An analysis of the numeri- 
cal results for the collider experiment STAR at RHIC shows that the radiative 
corrections prove to be important. 
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A Phase space of photon and scalar integrals 

The integral over phase space of the radiated photon can be presented in the 
form 

1 r rfik 

I[A] = - S[ip,+p,-h-knA] 

TT J ko 

where — i?„j/16 = Z\4 is the Gramm determinant, R^ = — A„jU^ — 2i3„jMi — C„j 
with the following coefficients: 

A„, = {s - vf - 2{s + v)mj + nil , 

Bui — [{s — f)(u — mg) + nil (2^ ^ t ^ f^i + 2to2)]w 

— [(s — nil )("^2 — t) + (u + nil — 2ni2 — 2?Ti; )v]z, 
Cux — [(v ~ ^ + fTi2)z + {u — nil )^] ~ ni2v[4:Z 

-2z(s -u~2v + m\- mf) + v(2u - mj + 2mf)]. (2) 

The limits of the double integration y™°^/™™ ^j^j ^max/mm ^^^.^ ^j-^^ roots of 
equations Ru^ = and u™°^ ~ u™*" respectively 



max /rain ^ui , / / -^Ui xn ^ui 

^vaaxlrmn ^ ]_,^ ^^ ^^2 _ 4^2^)^ T ^ S - V - m^ . 
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During the calculation we used the abbreviations 



T„ 



■m^,{v -u) + su, Tp^ = T^ + {v -u) m^r^, 



Tt - mlr{v-t) + st,T^r,^T^ + (v-tfmlrr^ 



w 



and following list of integrals over ui. 

_ 1 [ u\dui _ 3i?^^ - Au^Cui 



J2 = - 



where 



Js-- 



TT 



Jb-- 



Jr = - 



J7 = 



\/Rui 


2Al[^ 


uidui 


Bui 


\jRui 


.3/2' 


dui 


1 


\/Rui 


V ^«i 


dui 


1 


Ul y/Rui 


yC'ui 


dui 


Bui 


ufy/Ru, 


^3/2' 


dui 


1 


Z\^Rux 


yci:' 


du\ 


B., 



Bzi = -Au^ (v + z) - Bui , 



(3) 



(4) 



Cz,^AuAv + z)^ + 2Bu,{v + z) + Cu,. (5) 

At last, we list the following table of scalar integrals (VK-boson width is 
neglected where it is possible) 

h = /[I] = 1, 
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j^ ^ jX ^ log[(.s - vf/mfv] 



'2 l^2„ 
S — V 

2 /^2, 



r rr ^ 1 ^°S[(s + uy/mlv] 

h - J-[ — J ; , 

s + u 



Ul' 



^2 J 



U ^ mfl\^] 



s — V 



h = m^Ii^] 



1 , 1 



/e = I[ui] = -{s + u), 



.Ul uvil — log[(s — v)'^/mfv]) — st 

h = n — J = 7 ^2 ' 

z (s — vy 



h^i\ 



log- 



zzi vt m'lmf 



r 1 1 1 U 

h = I[ ] = log 2 2 ' 

zui vu m2mf 



Mt 



hn^Il^ 



■log 



(s-v)' 



rnfv 



z ' (s — vy 

s'^t^ + Auvst - Zv^v"^ 



2{s-vY 



hi = rnPQ] 



{s~ u)2 



(6) 



/i2 - Rel[ng] = - 



(s-v) 



log 



trv 



I,,^ReI[^]^^log 



^1 ^ n. 



i-W 



rp2 



mfmfytrwV 



/l4 = Rel[uinq] = -;r^ -^ H ^y- —3 log 



/i5 = Rel[ulng] = 



(s — w)^ (s — w)" 

2{s-vf 



trv 
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7 \A "^W 

2m'^twStvu — {sttw — rn^uvY 
{s — vY 



log 



ir-u 



"w 



he = Rel[—^] ^ hsimi -^ m2,u ^ t}, 

Ul 

In = rn^Reli-^] = , 



(7) 



hs = nngn+] = 



7r/2 — aictanit^ / {mw Tw)) ~ arctan(/V/mvF) 



mwrw{s — v) 



/lo = i\—n„n. 



Zl 






T„ 



Tfiwrw 



arctan 



{s - m'^)mw 



sFw 



arctan 



(s - m'^)Tu 



vsmy/rw 



+ {s + t) log 



mjTl 



mfm^v 



ho = I[ — n^n^] = /igJTOi ^ m2,t ^ m}, 



hi = mii[—{ngn+ - nin+)] = {s{v - 2m) + 2m^(« - 2v))snin+ /Tf^ 



h2 = mih^n„n+] = 



(s + u)^ 



(.2-tt 2^^9-'^9 
"1 



«T^* 



/23 = i[uin„n^^ 



1 






{st — uv) log — ^ 



"■w 



•■w 



Fw 



IT /2 — arctan 



l^L 



- 1 



hi = i[ulngn+ 



1 



(s — v)^ 



(vr — 2 arctan[t^/(mw/^M/)]) 
Amwrw 



X [3(sttu> — fn^uv)"^ — {st + m'^{s + u))^(s — u)^] 

/ i^ \ 
+s t { s — V ^- J — 4(s — w — t^)stuv ~ t^u V 
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+ [<u,s^t^ - 2stuv{t.u, + m^) + m^u^w^] log 

[v — u)^ 



Fw 



m 



w 






■log- 



+{vs-2{s + t)is-ml,))nin+\or "^ 

vs- {s + t){s -m^) 



ttiiyTw 



-is ^ ml,)nin+ 



( . As-'m'^)mw. . (s-mfy)Tu. 

X I arctani ■;^^ J — arctan[- '^ 



sFw 



he = I\ 






vsrav/Tw 



l2b{m\ -^ml.t^ m}, 



hi = I[ ] = — log- 2 2 ' 

ZlUl SV 171^1712 



i28 = I[ — \ = log 2 • 

Zl V — U TTlf W 



B Expressions for the Vj. 

B.l The case of qq — > l~Di,qq — »■ Z+i/j subprocesses 

Vi = Vi* = Auh -h- ul-! + ho - 2/ii 
Vig - Vil = Q,{-{2u{u -t)+ vt)l2 + iu- t){h - lu) 

+vt{2IIiu'^ - 2u + v)h - {t^{s + u) -vt~ 2u{u - t))/i2 
+ {t^{v -u){v- 2u) - t{2u^ - 2uv + v^))Ii3}/t^r 
+Qi'{2t^Ii + Au^h - 2t^ul3 - 2ulj - 2IIiu^vIg 
-2(t^ - M)((itu - 2u)Ii2 + hi - u{tw - u)Ii(i)}/twr 
-Qlu{vh^ + 2uhi} + Ql{-h + 2{nin+u^ - l)h + /i2 
+ (2t^ - m)/i6 + (itu - u)(4/i7 - t^/20 - 2{tw - u)l22) 
-t^hs} + QiQi'{{s - u)Iie - (s - u)/i8 
-(2m2 - 3uv + w^)/i9 - ((s - u){tw - 2u) + sv)l2o 

+ Sviv - 2u){l25 + h&) + 2sm2(/25 + I2Q - nin+ hv)} 



Va 



(8) 



(1) 



(2) 



(3) 
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Viw = {st^^Ii + 2sv?l2 + 2t-u]h - u{2s - v)I^ + (s - i;)/io 

-{a{tw - uY + u{tw{u -v) + su))Ii2 (4) 

-(tu. - u)(2itu + 2s - w)/i4 - (s + 2tw - v)Ii^}/twr 

Vqnj = Qt{sll2 - (2^2 - 3uV + V^)Il3 + (s(2u - tyj) + iP + vt))lis 
+ iitw -V + s)(2m2 - 3UV + V^) + SUv)Iig + t/23} 

+Qi,{-2h - {s - u)h + (2s + 4f^ - 4?i - v)h2 + 2/14 

+ {uv + 2{t^-u){s-u))hfi (5) 

+ {t^{Au + v -2t^- 2s) + u(2s - u - 3w))/i8 

— {twuv + (s - u)tw{tw — 2m) + 2u2(s — i;))/20 

+ (i - 2t^ + 2w)/23} 
K, = -s/i - /e + (2sf^ - u(3s + t))/i2 + (2s + 2t^ -u- v)hi 

+/15 + (i«,(u(3s + i) - st^) + u{uv -v^ - 2su))Iis (6) 

+ (u(2s - v) + tw{u + V - tyj - 2s))/23 - (s - i; + f^)/24 

B.2 The case of for qq — > l^vi,qq — > l~9i subprocesses 

Vi = v;(t^u) (7) 

^iq = ^g(* ^ u, mi ^ m2, Qi ^ Qi') (8) 

-2{t - t^fhl - /28} + Qlt{2tnin+h - Vho - 2t/22} 

+QiQi'{{s - t)/i3 - (s - v)hs 

(9) 
+ ((2t - t^)(s - t) - sw)/i9 - (2^2 - 3tv + w2)/20 

+ (2*2 - 2tv + v^)s{l25 + h& - nin+127) 

+nin+sv{v - 2t)i27} 

Vi^ = {2i„(i - v)h + {2st{t -v)- v^ujh + 2t^h 
+ {2st + tv + 2vu)Ij + (s - v)IiQ 

+ {{u - 3i„)w2 + 2{st + tl){v -t)+ U^{Av - i))/i2 (10) 

+ (2w(3i„ -u)-(t + t^){2t^ +v)- 2st)Ii4 
-{s -v + 2t^)Ii5}/t^r 
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V = ftK'S - 3i + 2w)/i2 + (2(s - t){t^ -t) + tv)Iu - 2/i4 
+ (f (2s - t - v) - t^is - 3t + 2v) + vu)Iis 
+ (t(2s + t^-v)(t^-t)-stl, + e{v-t^))Iig 

+ {2t^ -2t- U)l23 -{S- t)l28} + Qt'iiv - t)Ii2 

-(2f2 - 3tv + v^)Iie + (t(2s + t + t^,) + v{2u - i^))/i8 
+ {{s -v + t^){2e + v^) - tv{23 + 3(t^ - w)))/20 - u/23} 

K; = (t' - t)h -h + (2tt„(t - w) + t^ - 3to + 2v^)Ii2 

+ (t + 2t^-3v)Iu + Ii5 

(12) 
+ ((w - t - t^){v -t){v- t^) - s{2t{t -v)+ v'^))hs 

+ {{2S + t){v-t)-tl+ 3t^V - 2V^)I23 -is-V + t^)l2i 
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N,(P,) 




NjCPj) 




q(p2) 



piCk^r 



(1) 



(2) 



Fig. 1. (1) Sketch illustrating the interaction of parton 1 and 2 from incident 
nucleons Ni and N2, respectively. The partons carry fractional momenta xi and 
X2 and interact to produce a charged lepton and antineutrino . (2) The diagram 
for the lowest order subprocess qq ^ l^v. 




Fig. 2. Physically allowed region of variables xi and X2. Curves on the plot are 
X2 as a function of xi at different rj (kinematics of RHIC experiment STAR: 
VS=500 CeV, ki^^40 CeV). 
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Fig. 3. The virtual one-loop diagrams for qq -^ l^v process. The contributions 
to the self-energies and vertex corrections are symbolized by the empty loops, an 
explicit representation can be found in Ref.[17j. 







Fig. 4. Bremsstrahlung diagrams for qq ^ I i>"/ process. 
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dCT/dlk,-r (nb/GeV) 



0.05 




Fig. 5. Differential cross sections for pp -^ ji^X at \/S = 1.8 TeV, —1.2 < 
r] < 1.2, Z\^ = 27r (Tevatron) as a Junction o/fciy. Shown is the muon trans- 
verse momentum spectrum in the Born approximation (dashed line) and taking 
into account the total EWC (solid line). We use the MRS LO 98 set of parton 
distribution functions [21]. 
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;/dk,0/(d<To/dk,0 (%) 




27.5 



32.5 



47.5 



50 



37.5 40 42.5 

k.T- (GeV) 
Fig. 6. The ratio of the 0{a^) and the Born cross sections as a function of 
kiT at different values of quark masses: mu—5 MeV, md=8 MeV (curve 1), 
mu=30MeV, rud^SO MeV (curve 2), mu=100MeV, md=100 MeV (curve 3), 
mu=0.33 GeV , md — 0.33 GeV (curve 4-)- The rest of the parameters is identical 
to that in Fig. 5. 
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Fig. 7. Spin averaged part of the double differential cross sections for PP^ IX 
at \S = 500 GeV, A<I> — 27r (RHIC) as a function of ki^p for the pseudo- 
rapidity rj = —1 (curve 1), rj ^ (curve 2), rj — 1.5 (curve 4-), rj — 2 (curve 5) 
in the case of e , ji final states. Shown on the figures are the cross sections in 
the Born approximation (dashed lines) and taking into account the total EWC 
(solid lines). We use the GRV94 proton parametrization of parton distribution 
functions [22]. 
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Fig. 8. Corrections (5" /or PP^ l^X at RHIC kinematics as a function of kij' 
for the pseudo-rapidity 77 = —1 (curve 1), rj ^ (curve 2), rj = 1 (curve 3), 
rj — 2 (curve 5) in the case of e^ , fi^ final states. The rest of the parameters is 
identical to that in Fig. 7. 
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Fig. 9. Polarization 'part of the double differential cross sections for PP^ l^X 
at \/S = 500 GeV, A<I> = 2ti (RHIC) as a function of kij, for the pseudo- 
rapidity rj = —1 (curve 1), rj = (curve 2), rj = 1.5 (curve 4-), rj — 2 (curve 5) 
in the case of e , ji final states. Shown on the figures are the cross sections in 
the Born approximation (dashed lines) and taking into account the total EWC 
(solid lines). We use the GRSV96 LO proton parametrization of polarized parton 
distribution functions [23]. 
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Fig. 10. Single spin asymmetries for PP^ l^X at RHIC kinematics as a func- 
tion of kix for the pseudo-rapidity 77 = —1 (curve 1), rj = (curve 2), rj = 1 
(curve 3), rj ^ 2 (curve 5) in the Born approximation (dashed lines) and taking 
into account the total EWC: solid lines for the current quark masses, and dotted 
lines for the extreme choice rriu = rud = 0.33Gey. The rest of the parameters is 
identical to that in Fig. 9. 
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